The electron many-body problem in graphene 



Bruno Uchoa, James P. Reed, Yu Gan, Young II Joe, Eduardo Fradkin and Peter Abbamonte 
Department of Physics and Frederick Seitz Materials Research Laboratory, University of Illinois, Urbana, IL, 61801, USA 



o 

Oh' 



i 

s 

I 

o 
o 



> 
o 



X 



Diego Casa 

Advanced Photon Source, Argonne National Laboratory, Argonne, IL, 60439, USA 

(Dated: September 9, 2011) 

We give a brief summary of the current status of the electron many-body problem in graphene. We 
claim that graphene has intrinsic dielectric properties which should dress the interactions among 
the quasiparticles, and may explain why the observation of electron-electron renormalization ef- 
fects has been so elusive in the recent experiments. We argue that the strength of Coulomb 
interactions in graphene may be characterized by an effective fine structure constant given by 
Q*(k,cj) = 2.2/e(k, (jj), where £(k,tj) is the dynamical dielectric function. At long wavelengths, 
Q*(k,cj) appears to have its smallest value in the static regime, where a*(k ^ 0, 0) « 1/7 ac- 
cording to recent inelastic x-ray measurements, and the largest value in the optical limit, where 
Q*(0,a') ~ 2.6. We conclude that the strength of Coulomb interactions in graphene is not universal, 
but depends highly on the scale of the phenomenon of interest. We propose a prescription in order 
to reconcile different experiments. 

PACS numbers: 71.27.-|-a,73.20.Hb,75.30.Hx 



INTRODUCTION 

Graphene is a single atomic layer of graphite, that 
has been isolated only a few years ago|l|, and whose 
elementary quasiparticle excitations behave as massless 
Dirac fermions propagating in two spatial dimensions 
An important and to the present moment unsolved 
problem is determining to what extent electron-electron 
interactions are important in this material. In iso- 
lated free standing samples, a trivial estimate for the 
strength of electron-electron interactions in graphene can 
be achieved by computing ratio of the Coulomb energy, 
Ec = e^v}!'^ , where e is the electron charge, and n is the 
electronic density, to the Kinetic energy, Ek = hvu^l"^, 
where hv = 6 eVA is the velocity of the particles. Because 
the dispersion of the electrons is linear, in two dimensions 
those two energies scale in the same way with the density, 
and their ratio is a dimensionless constant known as the 
fine structure constant Q, 
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which is 300 times larger than the usual fine struc- 
ture constant in quantum electrodynamics, uqed = 

/{he) — 1/137, with c the speed of light. 

Although graphene is not as strongly interacting as 
other materials such as cuprates, where the strength of 
interactions is measured by the ratio U/t ^ 10, where 
U is the local Coulomb repulsion and t is the hopping 
energy, by this argument freestanding graphene should 
not be a weakly interacting system either. The validity of 
standard perturbation theory requires a < 1, which is not 
applicable in suspended samples. A coupling constant of 
K, 2.2 is large enough to result in a complete breakdown 
of perturbation theory. 



Why graphene fails to exhibit dramatic correlation ef- 
fects - even in two dimensions - is one of its most chal- 
lenging puzzles. Exacerbating the problem, a variety of 
experimental results that are sensitive to interaction ef- 
fects appear to give contradictory results, and further ap- 
pear to contradict what is expected from theory. Since 
the Fermi surface of neutral graphene (at half filling) is 
formed by points at the edges of the Brillouin zone, the 
Dirac points, metallic screening is not expected to influ- 
ence its electronic properties due to the vanishing density 
of states. At the same time, electron-electron interac- 
tions are generally expected to give rise to a logarithmic 
renormalization of the the Fermi velocity at the Dirac 
points 3, 



(2) 



where A is the bandwidth, which plays a role of an ul- 
traviolet cut-off, and q is the momentum measured away 
from the Dirac point. The structure of perturbation the- 
ory in graphene is such that this logarithmic divergence 
cannot be resummed in higher order expansion in a jj, 0] , 
in contrast to the behavior of Dirac fermions in one di- 
mension, the so called Luttinger liquids ,6]. In addition, 
the electron charge e does not renormalize, so the cou- 
pling constant a is expected to be logarithmically renor- 
malized to zero near the Dirac points, i.e. the interac- 
tion is marginally irrelevant. Although the damping of 
the quasiparticles is expected to be t{u)q) oc a'^\uJo\ to 
leading order in a, where |wo| — hvq is the energy of 
the quasiparticles, because of the logarithmic renormal- 
ization of V and a the ratio t(wo)/|wo| <SC 1 is small in 
the limit |wo| ^0, so the quasiparticles are still well de- 
fined (with logarithmic accuracy) near the Dirac points. 
In general, the renormalization of all physical quantities. 
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such as the compressibihty, susceptibihties, etc. can be 
derived directly from their scaUng dependence with the 
renormahzed velocity 

To be more concrete, one should expect for instance 
the single-particle spectrum of graphene, to be logarith- 
mically renormalized near the Dirac points. In supported 
samples, the fine structure constant is dressed by dielec- 
tric screening effects from the substrate, a = e'^/{hveo), 
with eo the dielectric constant. Angle resolved pho- 
toemission spectroscopy (ARPES) experiments do not 
show any evidence of a logarithmic renormalization in 
the spectrum (loj. even though in graphene on a mod- 
erate dielectric substrate such as Si02, where a ~ 1, a 
visible renormalization should be expected in the energy 
range of eV away from the Dirac point. 

Another interesting observable that can reveal renor- 
malization effects in graphene is the the electronic com- 
pressibility, K — (d'^E/dn^)^^ , with E the total energy, 
which measures the strength of interactions in an electron 
gas. This measurement has been recently carried out by 
a single electron transistor (SET) experiment 11| , which 
revealed large puddles of change in graphene at densities 
as low as lO^^cm"^, which corresponds to a Fermi level 
~50 meV away from the Dirac point. For non- interacting 
Dirac fermions, — v^J it /An. Many body effects were 
expected to give an additional logarithmic correction 
with the density due to the velocity renormalization 
The deviation from the non-interacting result, neverthe- 
less, was not observed. Although those experiments were 
not carried out on suspended samples, the modest dielec- 
tric screening from the Si02 substrate is not sufficient to 
explain the absence of verifiable renormalization effects. 
To make things more intriguing, recent ARPES measure- 
ments on doped graphene samples observed a non-trivial 
splitting of the Dirac cones [l^. which were attributed to 
the formation of plasmarons, a composite particle formed 
by an electron and a plasmon, the collective charge exci- 
tation of the Fermi sea. 

Many-body effects have also been observed in the pres- 
ence of strong magnetic fields [l3], which quench the ki- 
netic energy into discrete Landau levels. The most con- 
vincing evidence of interactions is based in the observa- 
tion of 1/ = 1/3 plateau in the fractional quantum hall 
effect in suspended samples 14|-|l6l| and in samples sup- 
ported on boron nitride 17[. At zero magnetic fields, on 
the other hand, there is no experimental evidence so far 
of the excitonic gap[l8j|. which has been predicted to open 
up in graphene when a > 1.1 [3]. In all current trans- 
port and spectroscopy exp eriments, graphene seems to 
behave as a semi- metal |2C 
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From the point of view of Coulomb scatterers, experi- 
ments involving the adsorption of K adatoms in graphene 
have reported a significant change in the electronic mo- 
bility of the samples [23|, what was interpreted as an in- 
dication that the transport in suspended samples was 
severely influenced by Coulomb scattering from charge 
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Figure 1: a) Energy spectrum of the electrons around the 
Dirac point, uj = hvq. b) Density map of the imaginary 
part of the polarization for non-interacting Dirac fermions, 
Imll'^' (g, cj). Dark blue color: Imn''^-'(q, cj) = 0; white re- 
gion: optical absorption edge at the border of the particle-hole 
continuum (uj > hvq). 



impurities. This interpretation is nevertheless at odds 
with another experiment 22] in which the dielectric con- 
stant of the Si02 substrate was enhanced up to two orders 
of magnitude by embedding the substrate in highly di- 
electric fluids. The variation of the mobility in this exper- 
iment was no more than 30%, indicating that Coulomb 
impurities do not influence the electronic properties of 
graphene. 

As we argue below, these apparent inconsistencies 
among different experiments and the difficulty of observ- 
ing electronic many body effects in the various physical 
observables may be due to the intrinsic dielectric screen- 
ing properties of graphene itself. In suspended samples, 
where a = 2.2, interactions among the electrons can be 
quite strong and lead to dynamical screening of the quasi- 
particles, which can be much more weakly interacting 
than previously believed. 



THE POLARIZATION OF THE VACUUM 

As the electrons interact, virtual processes that ex- 
cite electrons from the filled valence band up to the 
empty conduction band spontaneously create particle- 
hole pairs which eventually recombine and decay back 
into the ground state, ie. the vacuum. The process of 
spontaneous creation and annihilation of particles and 
holes polarizes the charge of the system. This polar- 
ization is expected to dress the quasiparticles and give 
rise to screening. Due to the absence of a Fermi surface 
(say at half-filling), the screening is dielectric in nature, 
in the sense that the Coulomb interaction remains long 
ranged but parametrically weaker. For a review about 
electron-electron interactions and the charge polarization 
of graphene, see Ref.|3|]. 

Although graphene is two dimensional, the fact that 
it is a semi-metal rather than gapped allows dielectric 
screening to emerge at long wavelengths. This is eas- 



ily understood if one computes the charge polarization 
function for non-interacting Dirac fermions in 2 + 1 di- 
mensions, 

U^'\q,L0)^-\ , (3) 

This polarization function is strictly real for uj < hvq and 
becomes purely imaginary for uj > Hvq, in the region of 
the particle-hole continuum where virtual excitation pro- 
cesses from the lower to the upper band are allowed (see 
Fig. 1). In the static regime {to = 0), the dielectric func- 
tion of graphene in random phase approximation (RPA) 
is 

e(g, 0) = 1 - F(q)n(i) (q, 0) = 1 -t- 7ra/2 « 4.45 (4) 
for a = 2.2^, where 

Vi,) = ^ (5) 

q 

is the Fourier transform of the Coulomb interaction 
V{r) — e^/r in 2D. For gapped graphene, as in any semi- 
conductor, there is a crossover in the behavior of the 
static polarization function at long wavelengths, hvq <^ 
A, n(^)((7, 0) — — g^/A, where A is the energy gap|25j. 
In two spatial dimensions, where the Coulomb interac- 
tion V{q) oc 1/q, one recovers the dielectric constant of 
the vacuum e(q A/hv,0) =^ 1, in contrast with 3D 
insulators iV{q) oc 1/q^), where the dielectric function 
saturates to a constant for q <C A/Hv. 

In the gapless case, if a = 2.2, there is no reason, a 
priori, to trust in the leading correction of perturbation 
theory in the physical observables, and one should seek 
to include diagrams of all orders in a. If one goes on 
to include the next correction to the dielectric constant, 
which is enclosed in the vertex correction of the bubble, 
V{q)n^^J,^^^{q,0) « -0.53a2. For a = 2.2, this term is 
almost of the same order of the RPA correction, in which 
case [2^ 

e(g, 0) = 1 + 7ra/2 + 0.530^ « 7 + 0{a^) . (6) 

In the opposite regime, where uj 3> hvq, one recovers the 
gapped situation, since the polarization bubble in leading 
order is 

nW(g^O,c.) = -i^, (7) 
4 lui 

and hence e{q — > 0,a;) — > 1. Therefore, the screened 
Coulomb interaction, and as a consequence the dressed 
fine structure constant of graphene, 

1 

a*(9>'^) = T — 7 (8) 

hv e(q, UJ) 
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Figure 2: Diagrammatic representation of a) the leading term 
in the polarization function, II'^' and b) the lowest order ver- 
tex correction. 

screening is ineffective and a*(0,a;) = 2.2 recovers the 
bare coupling constant of freestanding samples. In the 
on-shell intermediate case, where uj « hvq, the leading 
vertex correction to the polarization bubble has a loga- 
rithmic divergence near the on-shcU region (211]. 

Resummation of this divergence in the ladder channels 
to all orders in a has been proposed to give rise to a zero 
of the dielectric function corresponding to an excitonic 
bound state that lives in the optical gap of the particle- 
hole continuum [27,]. In any case, the RPA approximation 
seems to underestimate screening in the low energy sector 
UJ <^ hvq and does not account for excitonic effects at the 
edge of the particle-hole continuum. 

From the point of view of Coulomb impurities, the 
amount of charge induced by the polarization of the vac- 
uum around a test charge Q can be computed directly 
from the dielectric function assuming linear response, 

Qi„duccd('7) = -q(^1- , (10) 

where e{q,uj) = 1 — V{q)U{q,uj) is the dielectric func- 
tion. Since for Dirac fermions £(^,0) is a constant, the 
real space distribution of the induced charge is a delta 
function centered at the impurity. For a = 2.2, the over- 
all induced charge is Qj^ducod ~ —0.77(5 in RPA, and 
Qinduccd ~ —0.86Q + 0(q;'^) whcu thc first vertex correc- 
tion in the polarization is included. 



EXPERIMENTAL MEASUREMENT OF THE 
POLARIZATION 

One experiment that can provide direct information 
about the charge polarization properties of a system is 
the measurement of the response function with inelastic 
x-ray scattering. The response function is the imaginary 
part of the charge susceptibility, x(9, w), which is defined 
in terms of the full polarization function, Il(q,uj), 



has two distinct limits: the static regime, where 
a*{q,0) = (XG < 2.2, and the dynamic limit, where 



(11) 



4 



One major difficulty with this experiment for free 
standing graphene samples is that x-ray experiments can- 
not readily be performed on a one atom thick material. 
This difficulty can be overcome by realizing that although 
the susceptibility of the multilayer system can be quite 
different from the single layer case, they have very similar 
polarization functions at momentum and energy scales 
that are much larger than the electronic hopping energy 
between the different layers. The argument is the follow- 
ing: let us suppose for the moment that the interlayer 
hopping energy, is zero. The Coulomb interaction is 
long ranged and couples all the layers, 

with k = {kx,ky) the in-plane momentum. In the dia- 
gram in Fig. 2b for the vextex, it is clear that since the 
fermionic lines have no kz dependence (since the fermions 
do not disperse along that direction when t±_ = 0), the 
integration over k^ in the internal Coulomb line gives 

/.^oo 4^^3D(fc) = 27reV|k|, which is the Coulomb in- 
teraction of electrons in the single layer. This argument 
can be extended to all orders in perturbation theory in 
a. Hence, the polarization functions of the freestanding 
single layer and of the multilayer systems should be iden- 
tical provided they share the same single particle energy 
spectrum [soil . 

Since the external momentum and frequency enter as 
an infrared cut-off of the momentum integrals in the po- 
larization, provided that t±_ is much smaller than this 
energy scale, say t± <IC ma.yi{hvq,ui), t± can only give 
subleading corrections to the polarization but cannot af- 
fect the leading term. The similarity between the po- 
larization function of graphene and the multi-layer case 
breaks down if t± max(fiwg, w), in which case t± be- 
comes the infrared cut-off itself , giving rise to a crossover. 
The same rationale applies to other possible extrinsic in- 
frared energy scales, such as a Fermi pocket, Ep, which 
shifts the chemical potential away from the neutrality 
point, or an external magnetic field, B. Both are ex- 
pected to give rise to a crossover in the behavior of the 
dielectric function at sufficiently low energy scales, but 
do not affect the leading term of the dielectric function 
provided Ep, B <IC max(?iw(7, w). Hence, in spite of the 
fact that the dimensionality of the Coulomb interaction 
is different in the freestanding single layer and in the 
multilayer case, the polarization of the vacuum is quite 
similar in those two systems at energy scales where the 
single particle spectrum of the two is essentially the same. 

Examples of multilayer graphene systems include 
graphite in the Bernal AB stacking and other variations 
including turbostratic graphite, where the layers are ran- 
domly rotated. Information about the vacuum polariza- 
tion of the single freestanding layer can be also obtained 
experimentally from samples with a finite number of lay- 
ers. In single crystals of graphite, t± « 0.39eV is the 



interlayer hopping energy, below which the bands have 
hyperbolic dispersion near the K points of the Brillouin 
zone. In the turbostratic case, t± can be much smaller 
since the coherence of the hopping between different lay- 
ers is suppressed. Multilayer epitaxial grown graphene 
samples could also be suitable for x-ray experiments, and 
reveal detailed information about the vacuum polariza- 
tion properties of isolated graphene in the infrared. 

Recent inelastic x-ray experiments in single crystals 
of graphite conducted by the authors [111 have revealed 
that freestanding graphene is a highly polarizable system. 
Fig. 3 summarizes the result of those measurements. 
Fig. 3a shows the absolute value of the screened fine 
structure constant for a freestanding layer, which ranges 
from a*(q,0) ^ ac ~ 0.142 ± 0.092 « 1/7 in the static 
regime (uj <^ hvq) and long wavelengths, where screen- 
ing is more effective, up to uj) « 2.6 in the opposite 
limit, Lu ^ hvq, where screening is weak and collective 
modes such as plasmons emerge due to the diverging den- 
sity of states near the Brillouin zone boundary. The solid 
line LU = hvq is a guide to the eye, and indicates the edge 
of the particle hole continuum in the region of the Dirac 
cone {q < 0.5 A ). In panel 3c the static charge sus- 
ceptibility it is shown as a function of the momentum. 
This figure indicates that the leading behavior of x(9, 0) 
is linear in q, which is expected for a system with Dirac 
fermions and results in a finite screening strength at long 
distances. 

The imaginary part of the polarization function for the 
two lowest momentum data points, Imn((7,a;), is shown 
in Fig. 3b as a function of frequency and is compared 
with the imaginary part of the polarization for non- 
interacting Dirac fermions, Imn*^^)(g, w) (solid lines), de- 
fined by Eq. The peak in the sohd curves indicates 
the optical adsorption edge of the particle-hole contin- 
uum for a fixed momentum, depicted in Fig. lb. The 
data points show a redshift of the optical adsorption edge 
in ~ 0.6eV, which is similar in magnitude with the result 
of a recent ab initio calculation that explicitly accounted 
for interactions in the particle- hole channel[29|. This red- 
shift is interpreted as an excitonic shift in spectral weight 
due to Coulomb interaction between electrons and holes 
in the conduction and valence bands respectively. Hence, 
in spite of being a 2D semi-metal, graphene is polarizable 
to a degree similar to a conventional, 3D semiconductor 
such as Si or GaAs. 

Another way to view the significant, static polariz- 
ability of graphene is in the screening of a test charge. 
The reconstruction of the induced charge density in lin- 
ear response theory gives Qinduccd — — (0.924 ± 0.046)(5 
and the net impurity charge is only a small fraction of 
the charge of the original perturbation. In real space, 
the spatial distribution of the density is a local cloud 
of charge with a core of approximately Rq ^ 10 A of 
size, which quickly disappears beyond 15A away from 
the impurity (2S] . Since the Dirac Hamiltonian does not 
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Figure 3: (a) The magnitude of the effective, screened fine 
structure constant, ag{k,uj), reproduced from ref. [2^. (b) 
fmaginary part of the polarization function, H{q,Lo), mea- 
sured with inelastic x-ray scattering, compared with what 
is expected for ideal, Dirac Fermions (eq. 3). (c) Asymp- 
totic properties of the static, charge response function, x('?, 0), 
showing linear, leading behavior at small q. 



have a length scale, the corresponding charge response 
is purely local, as discussed below Eq. fTU^ . The high 
energy states by their turn affect the physics locally at 
length scales of the order the lattice spacing and cannot 
not directly influence the screening properties of the sys- 
tem at long distances. Therefore, in half-filled graphene 
the induced charge should be completely confined in a 
finite-sized region around the impurity. The implication 
is that the problem of screening a Coulomb impurity in 
graphene should manifest a complete separation of length 
scales. Any crossover induced by some infrared energy 



scale Eq, such as a small Fermi surface pocket, which 
restores metallic screening, or a small energy gap which 
leads to insulating behavior, will be manifested only be- 
yond length scales of the order of A ~ ttHv/Eo. For a 
small energy scale of Eq — O.leV measured away from 
the Dirac point, the corresponding length scale is of the 
order of A ~ 200A. Hence, in between these two different 
length scales, namely Rq and X ^ Rq, the behavior of 
graphene is quite universal and screens the test charge 
nearly completely. 

The value of the net charge found in the inelastic x- 
ray experiment suggests that graphene has an intrinsic 
static dielectric constant of the order of 10 [e{k,Q) w 15], 
which of course is just another way of visualizing the 
value given in the a; = 0, g — > region of Fig. 3a. This 
result implies that the transport is not significantly af- 
fected by Coulomb scatterers^2|. Although potassium 
atoms c hang e substantially the electronic mobility in 
graphene (23[, the scattering does not seem to be driven 
by the long range part of the Coulomb interaction, but 
rather by the short range part of it. 



PRESCRIPTION FOR PHYSICAL 
OBSERVABLES 

In suspended samples, the expansion in the unscreened 
a = 2.2 is poorly controlled,. The calculation of various 
physical observables such as the self-energy of the elec- 
trons and the total energy, from which one can extract 
the compressibility, can be instead organized in powers 
of the dressed Coulomb interaction, in the hope that 
the screening of the interactions among the quasiparti- 
cles will result in a convergent expansion. 

The use of the dynamically screened Coulomb interac- 
tion in the calculation of all physical observables may for 
instance reconcile apparent contradictions between dif- 
ferent kinds of experiments. For instance, the existence 
of plasmarons in graphene, which were proposed to exist 
in a region of the phase space where screening is typically 
weak (w 3> Hvq) , is not inconsistent with the observation 
that Coulomb impurities can be strongly screened out by 
the polarization of the vacuum, which is related to the 
static response of the system (w ^ fivq), where screen- 
ing is strong. The observation of the fractional quantum 
Hall effect is not inconsistent either with the idea that 
graphene is a strongly polarizable medium: the magnetic 
field in this case plays the role of an infrared cut-off that 
sets an energy scale below which the dielectric function 
undergoes a crossover in direction to a "gapped state". At 
energies comparable or smaller than this energy scale, the 
vacuum is strongly reconstructed by the Landau levels, 
which quench the kinetic energy and make interactions 
stronger. 

The lack of observable logarithmic renormalization ef- 
fects in the SET experiments on the compressibility in 
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the single layer supported on a Si02 substrate and also 
in ARPES experiments, can be attributed to dynamical 
screening effects, which inhibit the velocity renormaliza- 
tion at the energy scales currently reached by most exper- 
iments. Very recently, the measurement of Shubnikov de- 
Haas oscillations at very low magnetic fields {B ~ O.OIT) 
and densities as low as lO^cm"^ have revealed the first 
indirect observation of the logarithmic renormalization 
of the electronic spectrum in graphenejsij. An RG cal- 
culation incorporating dynamically generated screening 
within the RPA approximation yielded qualitative agree- 
ment with experiment. How to resolve this observation 
with x-ray experiments, in which the breakdown of RPA 
is very clear (Fig. 3b), is not immediately obvious. The 
assignment of an effective value for a in the Shubnikov- 
deHaas experiment {ac ~ 0.6 [29]) reflects an average 
over the strength of the Coulomb interactions in the en- 
tire dynamical range, and is not in principle inconsistent 
with the boundary values of a (1/7 in the static case 
and 2.6 in the optical regime) found in the x-ray experi- 
ment. It is perhaps suggestive that this effective value of 
aa ~ 0.6 approximately coincides with the value of the 
dressed fine structure constant at the edge of the particle- 
hole continuum [see Fig. 3 (a)], where on shell processes 
should dominate the renormalization of the spectrum. 
Further infrared experiments are needed to verify the ac- 
tual level of agreement between those two experiments. 
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Source supported by DEAC02-06CH11357. 



[1] K. S. Novoselov, et al., Electric field effect in atomically 
thin carbon films, Science 306, 666-669 (2004). 

[2] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, A. K. Geim, The electronic properties of 
graphene, Rev. Mod. Phys. 81, 109-162 (2009). 

[3] V. N. Kotov, B. Uchoa, V. M. Pereira, A. H. Castro Neto, 
F. Guinea, Electron-electron interactions in graphene: 
current status and perspectives, arXiv:1012.3484 (2010). 

[4] J. Gonzalez, F. Guinea, M. A. H. Vozmediano, Non- 
Fermi liquid behavior of electrons in the half- filled hon- 
eycomb lattice, Nucl. Phys. B 424, 595-618 (1994). 

[5] E. G. Mishchenko, Effect of electron-electron interactions 
on the conductivity of clean graphene, Phys. Rev. Lett. 
98, 216801 (2007). 

[6] T. Giamarchi, Quantum physics in one dimension, 
(Clarendom Press, Oxford, 2003). 

[7] 6 D. E. Sheehy, J. Schmalian, Quantum critical scaling 
in graphene, Phys. Rev. Lett. 99, 226803 (2007). 

[8] Y. Barlas, T. Pereg-Barnea, M. Polini, R. Asgari, A. 
H. MacDonald, Chirality and correlations in graphene, 
Phys. Rev. Lett. 98, 236601 (2007). 

[9] E. H. Hwang, B. Yu-Kuang Hu, S. Das Sarma, Den- 
sity dependent exchange contribution to dfi/dn and com- 



pressibility in graphene, Phys. Rev. Lett. 99, 226801 
(2007). 

[10] A. Bostwick, T. Ohta, T. Seyller, K. Horn, E. Roten- 
berg, Quasiparticle dynamics in graphene. Nature Phys. 
3, 36-40 (2007). 

[11] J. Martin, et al.. Observation of electron-hole puddles 
in graphene using a scanning single-electron transistor, 
Nature Phys. 4, 144-148 (2008). 

[12] A. Bostwick et al, Observation of Plasmarons in Quasi- 
Freestanding Doped Graphene, Science 328, 999 (2010). 

[13] J. C. Checkelsky, L. Li, N. P. Ong, Divergent resistance 
at the Dirac point in graphene: Evidence for a transition 
in a high magnetic field, Phys. Rev. B 79, 115434 (2009). 

[14] K. 1. Bolotin, F. Ghahari, M. D. Shulman, H. L. Stormer, 
P. Kim, Observation of the fractional quantum Hall effect 
in graphene. Nature 462, 196-199 (2009). 

[15] X. Du, L Skachko, F. Duerr, A. Luican, E. Y. Andrei, 
Fractional quantum Hall effect and insulating phase of 
Dirac electrons in graphene, Nature 462, 192-195 (2009). 

[16] F. Ghahari, Y. Zhao, P. Cadden-Zimansky, K. Bolotin, 
P. Kim, Measurement of the = 1/3 fractional quantum 
Hall energy gap in suspended graphene, Phys. Rev. Lett. 
105, 256805 (2010). 

[17] C.R. Dean, A.F. Young, P. Cadden-Zimansky, L. Wang, 
H. Ren, K. Watanabe, T. Taniguchi, P. Kim, J. Hone, 
K.L. Shepard, Multicomponent fractional quantum Hall 
effect in graphene. farXi^lOI071T79 (2010). 

[18] D. Khveshchenko, Ghost excitonic insulator transition in 
layered graphite, Phys. Rev. Lett. 87, p. 246802 (2001). 

[19] J. E. Drut, T. A. Lahde, Is graphene in vacuum an insu- 
lator?, Phys. Rev. Lett. 102, 026802 (2009). 

[20] Z. Q. Li, et al., Dirac charge dynamics in graphene by 
infrared spectroscopy. Nature Phys. 4, 532- 535 (2008). 

[21] 16 X. Du, I. Skachko, A. Barker, E. Y. Andrei, Approach- 
ing ballistic transport in suspended graphene. Nature 
Nanotechnol. 3, 491-495 (2008). 

[22] L. A. Ponomarenko, et al, Effect of a High-K Environ- 
ment on Charge Carrier Mobility in Graphene, Phys. 
Rev. Lett. 102, 206603 (2009). 

[23] J. H. Chen, C. Jang, M. S. Fuhrer, E. D. Williams, M. 
Ishigami, Charged Impurity Scattering in Graphene, Na- 
ture Physics 4, 377 (2008). 

[24] Of course, since a runs to zero logarithmically, one should 
eventually expect e{q, 0) ^ 1 at very long wavelengths. 

[25] V. N. Kotov, V. M. Pereira, B. Uchoa, Polarization 
Charge Distribution in Gapped Graphene, Phys. Rev. 
B 78, 075433 (2008). 

[26] V. N. Kotov, B. Uchoa, A. H. Castro Neto, Electron- 
Electron Interactions in the Vacuum Polarization of 
Graphene, Phys. Rev. B 78, 035119 (2008). 

[27] S. Gangadharaiah, A. M. Farid, E. G. Mishchenko, 
Charge response function and a novel plasmon mode in 
graphene, Phys. Rev. Lett. 100, 166802 (2008). 

[28] J. P. Reed, B. Uchoa, Y. II Joe, Y. Can, D. Casa, E. Frad- 
kin, P. Abbamonte, The effective fine structure constant 
of freestanding graphene measured in graphite. Science 
330, 805 (2010). 

[29] 25 L. Yang, J. Deslippe, C.-H. Park, M. L. Cohen, S. 
G. Louie, Excitonic Effects on the Optical Response of 
Graphene and Bilayer Graphene, Phys. Rev. Lett. 103, 
186802 (2009). 

[30] Although the argument considered a continuum of lay- 
ers, it can be easily generalized to the discrete case as 
well, by noticing that the Coulomb interaction in this 



7 



situation is given by Vsnik) — 27re^dS'(fc)/|k[, which 
is the Fourier transform of the e^/r interaction for an 
infinite stack of layers separated by a distance d, and 
S{k) = sinh(|k|d)/[cosh(|k|d) — cos{kzd)] is the structure 
factor. Integration over in the Brillouin zone of size 



27r/d, namely /^^^^ ^Vaoik) = 27re^/|k|, which again 
is the Coulomb interaction in 2D. 
[31] D. C. Elias et ffl;.. larXiv:1104.1396r 2011). 



